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TECHNICAL NOTE 2869 

REFLECTION OF WEAK SHOCK WAVE FROM A 
BOUNDARY LAYER ALONG A FLAT PLATE. II - INTERACTION 
OF OBLIQUE SHOCK WAVE WITH A LAMINAR BOUNDARY LAYER 
ANALYZED BY DIFFERENTIAL-EQUATION METHOD 
By Yung-Huai Kuo 

SUMMARY 


By analogy with the boundary-layer concept, the flow produced by 
the interaction between a shock wave and a laminar boundary layer is 
subdivided into a viscous layer and a potential field. The assumptions 
that the compressibility effect in the inner layer is negligible and 
that the original flow in the outer layer is uniform lead to simple 
analytic solutions for the perturbed flow. The joining conditions at 
the interface between the layers determine an eigenvalue which gives 
the rate of decay and the character of the disturbances both upstream 
and downstream of the point of incidence. The final conclusions are 
in agreement with experiments. 


INTRODUCTION 


The present investigation is an independent study of the inter- 
action of an oblique shock with a laminar boundary layer in a compress- 
ible supersonic stream. In reference 1, where interaction of weak shock 
waves with both laminar and turbulent boundary layers was treated, the 
integrated momentum across the boundary layer was considered, rather 
than the balance among various dynamic effects at each point. This 
momentum-integral method is simple and, in certain respects, powerful 
and capable of yielding useful qualitative information such as the 
upstream pressure influence, pressure distribution, and the gbowth of 
boundary-layer thickness due to the presence of a shock, but it fails 
in regard to what actually happens inside the boundary layer. In the 
present report a different approach has been adopted, with the inten- 
tion of filling the gap left by the previous investigation. The pur- 
pose will, on the whole, be complementary, so as to provide a physical 
picture for the understanding of this complex phenomenon. 
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Contrary to reference 1, the differential-equation method is 
employed here. According to available experimental observation, when 
an oblique shock is incident upon a laminar boundary layer the result- 
ant flow bears no resemblance to the flow predicted by potential theory. 
For if the viscous flow is absent the flow ahead of the shock will not 
be affected. Because of the presence of the boundary layer in which 
there is a subsonic layer, however, a sudden decrease of pressure at a 
point will immediately be transmitted forward by the inability of the 
subsonic layer to support an excess pressure rise. When the pressure 
is transmitted, the flow in the boundary layer will be retarded and the 
streamlines distorted. Since the outer field is supersonic, this change 
occurring in the viscous layer will affect the whole potential field. 

This is actually observed. For stronger shocks, the flow in the bound- 
ary layer generally will separate and will have backflow under the 
influence of an adverse pressure gradient. An adequate theory that is 
able to account for the observed effects cannot be formulated unless 
the boundary-layer hypothesis is abandoned entirely. One therefore is 
faced by a much more difficult mathematical problem. 

To restrict the scope and complexity of this study, let it be 
assumed first of all that the boundary is an insulated flat plate, and, 
secondly, that the incident shock is weak and its angle of incidence is 
such that regular reflection would be possible, were the flow friction- 
less, and, lastly, that the free-stream Mach number is not large. Under 
these assumptions, the angle of deflection of the flow in passing through 
the shock wave will be small, and the temperature variation between the 
free-stream condition and that of the plate will not be large. In fact, 
the study of a laminar boundary layer indicates that, for moderate Mach 
number, the temperature as well as the compressibility effects are unim- 
portant (reference 2). This must remain true even if the flow is not 
boundary-layer flow. Therefore, without loss of generality, it will be 
assumed that the viscosity and thermal conductivity of the gas will be 
taken as constant and the Prandtl number is unity. 

In order to bring the interaction problem within the scope of prac- 
tical mathematical analysis, these simplifying hypotheses have to be 
made in the absence of a proper method of approximation, such as the 
boundary-layer theory. Broadly speaking, examination of schlieren 
photographs of the flow produced by the interaction between a shock and 
a laminar boundary layer will reveal that two characteristically dif- 
ferent outer and inner regions exist for sufficiently high Reynolds 
number. The outer field is characterized by its strong potential char- 
acter, whereas the region close to the wall is predominantly viscous, 
which is quite reminiscent of the boundary-layer flow. It appears 
natural, therefore, to assume a priori that the viscous effect is con- 
fined to a thin layer in the vicinity of the boundary and the outer main 
flow is potential. These two different flows are then in dynamical 
equilibrium. If one is disturbed, the other will be affected. Since 
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the outer field is supersonic, any local change will he felt in a much 
larger region than that in a subsonic field. 

After the flow field is separated into two regions, specific assump- 
tions regarding the structure of the viscous layer can be made. It is 
important to note that in the case of a compression incident wave, the 
overwhelming effect taking place in the viscous layer is the sudden 
decrease of the velocity or even reversal of flow (backflow). If back- 
flow sets in, the flow speed in the subsonic region will be very much 
reduced. As a result of this, the streamlines will be pushed outward 
and the flow compressed. On account of the displacement of streamlines, 
the subsonic region will become thicker, and the thickening of the sub- 
sonic region is a decisive factor that distinguishes the strong from 
weak interactions. The importance of this dimension (the thickness of 
the subsonic layer) has already been established by Tsien and Finston 
(reference 3 ) in the case of inviscid theory. 

The viscous layer, thus, must have two distinct sublayers, each 
displaying a different character. In the supersonic layer, the flow is 
characterized by large velocities and is not unlike that of an unsepa- 
rated boundary^ layer. Therefore, for this layer both viscous and inertia 
forces are important. On the other hand, in the subsonic layer, espe- 
cially with backflow, the average speed will be very small. In this 
case, because of slow motion and moderate temperature change, the change 
of density is always a lower -order effect. In fact, the contribution 
due to compressibility is proportional to the square of Mach number, 
and, if the average Mach number in the subsonic region is small, the 
compressibility effect is, indeed, negligible. Because of this physical 
fact, the subsonic layer will be taken as incompressible. 

With these assumptions, the problem is finally solved by pertur- 
bation of weak incident waves. As a test of these assumptions, a simple 
flow with broken-line velocity profile is taken as the basic flow: In 

the incompressible layer, the velocity is a linear function of the dis- 
tance from the plate} in the compressible layer, it is constant. The 
density in the basic flow is constant in each layer but takes different 
values. Thus, at the interface where the two flows join, the velocity 
is continuous but density is discontinuous. For this case, a first-order 
solution consistent with these assumptions is completely determined. 

In the case of weak shock, excellent agreement with experiments has 
been achieved for the pressure distribution on the plate. It confirms 
the conjecture that separation of the flow as well as backflow always 
occur. Because of the occurrence of backflow, transition downstream of 
the point of incidence is unavoidable in the viscous layer. There are 
strong experimental evidences but detailed investigations are yet to 
be conducted. In the outer field, on the other hand, it is predicted 
that in the place of the regularly reflected shock there is a strong 
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expansion, and farther downstream a train of strong compression waves 
must exist, eventually forming an envelope. Therefore, downstream of 
the point of incidence a second shock must occur. This deduction is 
also confirmed by experiments. 

Lastly, the importance of nonlinear effect is discussed. 

This study was conducted at the Graduate School of Aeronautical 
Engineering, Cornell University, under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
The author wishes to acknowledge the efficient assistance of Mr. Nelson 
H. Kemp. 


SYMBOLS 


A constant 

a' speed of sound 

C 1 ,C 2 ,C 2 ,C^,C 5 ,Cg constants 

D 
d 

d/dt 


constant 

constant 

convective derivative 


F,G 


f l'«l 


H 



J l/3 

k 1 ,k 2 ,k 3 ,k 4 


scalar functions 

defined by equations (43) and (kl), respectively 
nondimens ional enthalpy 

Bessel function of first kind with imaginary argument 

Bessel function of first kind with real argument 

» 

constants 
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Bessel function of second kind with imaginary argument 


k = (7 - l^/R 2 


k = 

L 

M 

n 

P 


P l> P 2> P ] 

P 

P' 


1 + i( 7 .i)?^ 
3 R 


1 + 4 7M 2 xA 
3 R / 


R 


length of plate from leading edge to point of incidence 

Mach number 

dilatation 

pressure 

constants 


nondimens ional pressure 
pressure 

Reynolds number (u^L/v^) 


(p’/poo^ 2 ) 

1/2 


= r(i/3)/: 


>2/3 


T 

T' 

s,t 

U 

U, V 
U ' , V ' 

u (l) . 


nondimens ional temperature 
temperature 

defined by equations (B4) 
velocity 

nondimens ional velocity components 

velocity components 

constants 
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X - -a/xg 
x,y 
x',y' 
z(n) 

z 

a = O.332 
a = P 0 (°W 

P = \/m 2 - 1 
Pi 

P 2 


nondimens ional Cartesian coordinates 
coordinates 

defined "by equation (CA) 
defined by equation ( 30 ) 


defined by equation (A 3 ) 
defined by equation (C15) 


r (z) 

7 

A 

5-3 

e 

T) 

n 

A (r) ) 

\ 

V 

Y 00 

P 

P* 


gamma function 

ratio of specific heats 

Laplacian operator 

deflection of flow (equation (B 6)) 

flow-deflection angle 

defined by equation ( 31 ) 

defined by equation (C6) 

defined by equation (C5) 

eigenvalue 

free-stream kinematic viscosity 
nondimens ional density 
density 



MCA TN 2869 


7 


P Q density at y = 0 

a = XR - X 2 

t = \(\a) ^ ^ 

r defined by equation (C7) 

cp nondimens ional velocity potential 

\|r stream function 


0) 

Subscripts: 

C 

c 

i 

0 

P 

s 

sep 

t 

x 

y . 

1 

2 

3 


angle between velocity and shock 
complementary 

solution in compressible layer 

solution in incompressible layer 

value on plate 

particular 

due to step wave 

at separation point 

due to transmitted wave 
partial derivative with respect to x 
partial derivative with respect to y 
in region 1, before shock (see fig. l) 
in region 2, behind shock (see fig. l) 
in region 3 (see fig. l) 


00 


free stream 



8 


MCA TN 2869 
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parallel to oblique shock wave 
perpendicular to oblique shock wave 


Superscripts: 

(0) 

(1) 

(2) 

(3) 


zeroth order 
first order 
second order 
third order 


STATEMENT OF PROBLEM AND BASIC ASSUMPTIONS 


Let there be a laminar boundary layer in a compressible viscous 
fluid along an insulated flat plate immersed in a steady uniform super- 
sonic stream, and let an oblique shock be incident upon the plate* When 
the steady condition is established, the flow in the neighborhood exhibits 
a character which is entirely different from the original flow* 

From all indications, this flow does not obey the boundary -layer 
approximation) nevertheless, for simplicity for future analysis, the 
concept of boundary layer, or viscous layer, will be retained* Namely, 
the whole flow field is visualized as consisting of inviscid and viscous 
flow in equilibrium with each other. The possibility of existence of 
such a demarcation line will be assumed at this moment* The solutions 
obtained are consistent with the assumption, as will be seen later, so 
that the theory is self-consistent at least* Naturally, its further 
justification rests upon experimental evidence. 

The main feature in the viscous layer is that, in the case of a 
compression wave in the inviscid outer flow, backflow generally exists. 

For this reason, terms in the equations of motion which are unimportant 
according to boundary-layer approximation become decisive as the sup- 
posed large-order terms vanish. Therefore, the pressure gradients along 
both directions have to be considered. To simplify the mathematical 
process, some minor effects, such as the variation of the viscosity 
coefficient and thermal conductivity with temperature, will be neg- 
lected and the Prandtl number will be taken equal to 1. For moderate 
Mach numbers, less than 3> say, this neglect, according to boundary- 
layer investigations, will have little effect on the major character- 
istics of the flow* 

The inviscid flow generally is rotational, as it involves shocks. 

This is particularly true in the case of a local supersonic zone over 
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a curved surface. For a purely supersonic flow, if the shock is 
slightly curved, the vorticity generated behind the shock is of high 
order and can be neglected. The perturbed flow in the outer region can 
then be regarded as irrotational. 

As the region which is influenced by the presence of the shock, 
according to experimental observations, is confined to an area about 
the point of incidence, with a dimension only a fraction of the length 
of the plate, a point which lies at a distance about the length of the 
plate from the point of incidence will be considered as at infinity. 
Consequently, the boundary-layer flow will be replaced by a "shear flow" 
extending to both positive and negative infinities. This approx ima tion 
is justified if the derivative along the flow is much larger in the 
perturbed flow than that in the original flow. For large Reynolds num- 
ber, this condition can always be satisfied. 


METHOD OF SOLUTION 


The flow is supposed to be two dimensional and steady and the fluid 
is compressible. If the flow field can be subdivided into viscous and 
nonviscous regions, the flow in the viscous layer satisfies in dimension- 
less variables the system: 

0 5? ' + i( Au + 1 “*) (1) 

/ 

p S “ - p y + §( Av + 5 ”y) < 2 > 

(PU) X + (pv)y =0 (3) 


p 22 = i ah + lr~:- 1)M2 (^) 

dt R R \dt/ x 





( 4 ) 


and in the case of perfect gas 
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Here the subscripts denote partial differentiations with respect to the 
Cartesian coordinate x or y indicated; the Laplacian operator A , 
the convective derivative d/dt, the nondimens ional enthalpy H, and 
the dilatation n are defined by 



Furthermore, the velocity components u' and v', the pressure p', 
the density p', the temperature T', and the coordinates x’ and y' 
are nondimens ionalized as follows: 
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where p m , and v M are free-stream velocity, density, and kinematic 

viscosity, and L is the length of the plate from the leading edge to 

/ x 1 / 2 

= (Hz*) 


the point of incidence. Finally, the parameters 7, M, and R 


XVoo / 


stand for, respectively, the ratio of specific heats, the free-stream 
Mach number, and the Reynolds number. 


On the other hand, the inviscid flow in the outer region is assumed 
to he irrotational and, consequently, is determined by a nondimens ional 
velocity potential cp(x,y ) satisfying the equation 


(a 


) 2 - 


(u') 


cp. 


XX 


2u*v'cp. 


xy 


(a'f - (V)‘ 


cp. 


yy 


= 0 


(8) 


where the sound speed a' is given in terms of the velocity components 
u' and v' by the relation 

(a'f - (a„' f - lUl £»• f + (y f - u/j (9) 


By the assumption that these two different flows are in equilibrium with 
each other, the flow determined from equation ( 9 ) must join smoothly with 
that given by the system of equations (l) to (5) subject to the boundary 
conditions 


when y = 0 and 


u = v = 0 


u = cp x 

V = CPy 


when y = 00. 

As the disturbances are initiated by the incident shock, it is 
expected that they are small if the incident wave is weak. From experi- 
ence, this is at least the case in the field upstream of the point of 
incidence. If the shock strength is characterized by the flow-deflection 
angle e, then for small values of e the solutions are expanded in 
powers of e : 
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u = + su(^) + e 2 u( 2 ) + 

v = + g2 v (2) + 

p = p(0) + e p(l) + e 2p(2) + 

p = + ep^^ + e^p(2) + 

T = + e T^ + e 2 T^ 2 ^ + 




( 10 ) 


By substituting the quantities defined in equations (10 ) into the system 
of equations (l) to (5), there results the following system of equations, 
according to the powers of ej namely, for the zero order. 


(0) a ( 0 ) u (0) _ (0) 1 

dt x R 


(o) a<°>v<°> 


at 


- - P <°> + i 

y r 


|Au<°) + i ^( 0 ) 


. ( 0 ) 1 ( 0 ) 

Av + O “y 


(p(O)u(O)} + (p< 0 ) v(°>) .0 
x y 

(0) a(°>H (0) _ 1 (0) + (7 - Dm 2 

dt R R 


> ( 11 ) 


U< 0 V°A + (d< 0 >v<°> 


dt 


dt 


£(»„«») - 2(v<0) n (0)) 

3 x 3 y 


p(°) = ( 7 M 2 )~ 1 p(°^T(°^ 
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where 


d^ (0) 8 ( 0 ) 8 

dt 8x 8y 

H (°) = T^ 0) + 2LzJ; M 2 ^ 0 ) 2 + v (°) 2 ) 


> 


- (0) - ^ ♦ v y C) 


for the first order ? 

p(0)Wil + W + D (1) d(°) u <°) . 

' dt dt / dt 

|( Au(l) + ^ % (1) ) 

p (0)(dlS^ + d(l) v (0)> ) + D (l) d<°V°> 

\ dt dt / P dt 


-P (1) + 


X 


-p « ♦ 


1(av (1> + 1 n y V) 


(p(O)u(D + p(Du(0)) + (p(0) v (l) + p (l) v (0)) , 

X ' v 


' (0) ♦ (p (1) ^ - p (0) = ± *,(i> + 

a t \ dt dt / R 


(7 - l)M^ 


R 


+ d^u^' 


lO-KW 


dt 


dt 


dt 


dt 


|(u(0)n(D + u(Dn(0)) - |(v(°) n a) + T (l)„(0)) 
^ x 3 } 


p(l) = (^M 2 ) + p( d )i(0)) 


J 


( 12 ) 


> ( 13 ) 
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where 


,(D 


d — - A 1 * A + v d) A 

dt & ^ 


H 


(1) = T (l) + (7 - 1) u (0) u (l) + v (0) v (l) 


> 


n(D = % (D + v (D 

«/ 


and for the second order ^ 


(14) 


.Co) 


iAA + AA (1) + dgjs ft>) + p(1 )(d^«w , d< l V0)\ d (0) u ( 

' dt at it r p \ dt + at J + p ST 


,(o) 


- p x (2) + |( Au(2) + J >\ (2) ) 


.r 0 )(d (0) v (2) d^V 1 ) 

P \ r: + rr + / + p 


dLt 


dt 


dt 


(i)/d (0 V (1) + d^V^j + ( 2 ) a 

V dt dt / 


(0) T (0) 


dt 


(p(0)u(2) + p (l) u (l) + p (2)u(0)^ * ( p (0) v (2) + p(l) v (l) + p (E) v (0)J , Q 

p ( o, l^ll . (pC) ill . .« ig^W * ( B (o) ig . p(» 11 ; 


,( 0 )’ 


> (is) 


^)h(o) 

dt 


J AH^ + 

1 

H 

K 

/d(°)u( 2 ) . d 

R 

R j 

[V dt 


dW 1 ) 

+ d ( 2 M°)\ 

\ dt 

+ dt 

dt Jy 

|( v (0) n (2) 

+ tW 1 ) + 


y 3 


P (2) = (rM 2 )‘ 1 (p ( 0) T (2) + + p( 2 V 0 )) 
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where 

d(2) 

at 


u(2) A + t (2) A 

dx Sy 




h (2) = t ( 2) + ( 7 _ i)j^^i(°) u ( 2 ) + v( 0) y( 2 ) + — ±_zliL^ 


n< 2 ) - ^(2) + v„(2) 


y 


In the case of the potential flow, if the velocity potential is 
expanded as 


cp = x + <sp 


( 1 ) 


e V 2) 


( 17 ) 


there will occur, similarly, equations for the first- and second-order 
quantities. These are: For the first order 



( 1 ) 


<P- 


( 1 ) _ 


yy 


= 0 


( 18 ) 


and for the second order 


P 2 cp. 


XX 


(2) . 9 .J 2) . _ 2m 2 


yy 


(! ♦ H 1 ^K ( 1 ) <Pxx a) - V 1 V 1J 


( 19 ) 


where p 



FIRST-ORDER SOLUTION - UPSTREAM, x < 0 


Let the point of incidence he chosen as the origin of the Cartesian 
coordinates. Then, negative x will he called upstream, and positive x, 
downstream, of the point of incidence. The various regions will he num- 
bered hy 1, 2, and 3 as shown in figure 1 in which OS indicates shock 
and OM, the limiting Mach line of region 2. 
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As it has been assumed that the basic flow is a laminar boundary- 
layer flow of a Blasius type; such as considered by Von Karman and Tsien 
(referenced) and by Emmons and Brainerd (references h and 5); system (ll) 
simply reduces to the well-known Prandtl boundary-layer equations; whose 
solution has already been found in these references. 

Since the basic flow is a boundary -layer flow with constant pres™ 

sure system (13) can be simplified* From experiments it has been 

established that the space variation of the perturbed quantities is much 
more rapid than that of the unperturbed quantities. By the boundary- 
layer approximation; the ratio of the partial derivatives 8/cht: d/dy 
is 1:R. For large Reynolds number; the x-derivatives of the basic 

flow; as well as v^, which is of the order R“^; can be neglected in 

comparison with the y-derivatives* That is, u^ and are func- 

tions of y only. This approximation is confirmed by experiments and; 
as a matter of fact; it is customarily used in the pressure measurements ; 
because the wave is shifted forward and back relative to a fixed pressure 
orifice to measure the pressure distribution before and behind the wave. 

In the case of unit Prandtl number; is a constant. Then; by 

dropping terms such as v^^Uy^) and u x ^^u^^, system ( 13 ) becomes 


p (0)^ u (0)u x (l) + Uy(0) v (l)^ ss -p x ^ 1 ) + i^Au^) + i n^ 1 ^ 
p(°) U ( 0 )v x ( 1 ) = 4 * ^Av^ 1 ^ + ~ ny^ 1 )^ 

p( 0 )n(-^) 4 - u^^p 4 - v(^-)p = 0 

x y 



(20) 


= (yM 2 ) 1 (p^°' ) T^ 1 ^ + p^V^) 



NACA'TN 2869 


IT 


Incompressible Layer 

It is noted that the coefficients of system (20) depend on the 

basic velocity profile u/°^(y). Since u^(y) cannot be expressed 
by simple functions, in order to simplify the analytical work further 
simplifications are necessary. 

First of all, in the case of insulated plate, the temperature and, 
hence, the density in the basic flow will have a vanishing gradient on 
the plate. This makes the variation of the density in the viscous layer 
much smaller than that of the velocity. When the flow is subject to an 
adverse pressure gradient, the flow will be further retarded. When the 
backflow occurs , the average speed will be very low and, consequently, 
the representative average Mach number will be small. Under this condi- 
tion, the change of density is less important than that of the pressure. 

A layer where this approximation is valid is called an "incompressible 
layer" and the dominant effects will be pressure and frictional forces. 

Of course, it is difficult to define the thickness of this layer 
beforehand. Generally, it would correspond to the subsonic portion of 
the viscous layer. In the basic flow, the sonic boundary can be exactly 
calculated. When the flow is perturbed, it is unknown, but certainly 
will be thicker, for the flow is subject to an adverse pressure gra- 
dient. Owing to this fact and also for mathematical expediency, the 
thickness will be taken in accordance with the way the velocity profile 

u (°)(y) is represented. This will become clear below and the thickness 
defined turns out to be nearly half the original boundary -layer thickness. 

Because of the assumption that density is constant and assumes the 
value on the plate, system (20 ) simplifies to 


- o (0) ( u(0 V 1) + U v ( 0,t(1) ) * -Px (l) 




p (o) u (o) ( 1 ) , ( 1 ) + 1 AV (1) 

O -X. J R 

nd) = 0 




(o) u (o) H ( 1 ) .land) + ? (r i)m2 uJ°)v d) 

> x r r y 


X 


(21) 
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Thus, the velocities, pressure, and temperature can he dealt with inde- 
pendently. For the velocities, the first three equations yield by 

elimination of p 



After the velocities are known, the temperature will be given by the 
energy equation. In the present problem, however, temperature is not 
an interesting quantity and will subsequently not be mentioned. To 

solve equations (22), let a stream function i|f^^(x,y) be introduced 
.by the relations 

ud> » v (1) 

v^ = -\|r x ^ 

From equations (22) i|r^^ then satisfies the equation 

i A Ai|r ^ ^ + P o ^°^Uyy(° (23) 

By neglecting the compressibility effect in the inner layer, the 

basic velocity profile u^°^(y) becomes the Blasius profile. The vari- 
able coefficients of equation ( 23 ), for small values of y, will be 
power series in y. However, it has been recognized that in the case 
of Blasius profile the velocity attains the free-stream value fairly 
rapidly and the initial portion is nearly linear. Consequently, it will 
not involve serious error to replace the continuous profile by a broken- 
line one, so that the initial part is proportional to y and the 
remaining part, constant with the free-stream value unity. If the skin 
friction agrees with the exact value, the velocity profile can then be 
defined as 
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If the velocity is continuous at y = y^, then y = a" 1 which is about 
half the houndary-layer thickness* namely y ~ 3- If the incompressible 
layer is defined as the interval 0 = y = y-, * then equation (23) becomes 


ay A|r ^ = A 


(25) 


where a = p ^°^aR. 
o 

To solve equation (25 ^ it is natural to assume the form 


^( 1 ) = t(y)e Xx 


X > 0 


1 


( 26 ) 


for region 1 where x < 0. This assumes a special form of compression 
waves induced in the potential field; that is, the solution of equa- 
tion (l8) 


cp (l) = A e X(x ' Py) (27) 

X 

where A is a constant to be determined. This is entirely in agreement 
with the conclusions reached in reference 1. Substituting equations (26) 
in equation (25 ) and simultaneously putting 

t + X. 2 t = z (y ) ( 28 ) 

there results 

z 1 * + (lay - l 2 )z = 0 (29) 


The general solution of this equation is 
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where 

ri = (Xa) 1 ^ 3 (y - y&T 1 ) ( 31 ) 

C-l and Cg are constants, and 1 ^ 2 . ^ 3 / 2 ^ k^S, t)^/ 2 ^ denote the 

3 3 

Bessel functions of the first and second kind with imaginary argument. 
With z known, the solution of equation ( 28 ) is 

t = Cj cos tt) + sin ttj + 

8111 — f (cos ttj)z(t}) dr) - — f (sin tt))z(tj) d-q ( 28 a) 
‘ J % T J *lo 


with t = l(Xxi) - ' 1 ’/ 3 - On the plate y = 0 and = v^^ = 0; hence 

t(0) = ^'(0) = 0. For large Reynolds number and small values of X, 
y = 0 corresponds approximately to T] = 0. The boundary conditions 
thus require = 0. Moreover, at y = y^, rj will he large, 


and since ttj^ = Xy-|_< 1, the integrals J' -j^ . TqJ>I^ ^ r j^/ 2 ^ dq 

pr\ | q3/2 

will diverge like / dq as q approaches infinity. Therefore, 

Jo 

C% = 0. The solution in the incompressible layer can then be read as 
follows 
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Thereby the velocity components are 



On the other hand , the pressure can be obtained by integration of 
the first of equations (21 ) and is, by the form of solution, 

p(D = _p <o) u (o) u (D . Po (o y°h (1 > jA( 1 ) 2 

o \ XR X 7 

where the constant of integration vanishes by the condition at negative 

infinity. Substituting u^ and v^ 1 ^ from equation ( 33 ) and making 
use of equation (28), a straightforward reduction yields 
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On the plate y = tj = 0 this reduces to 


P (l) (x,0) 



(3^a) 


Compressible Layer 


In the compressible layer, because of the forward momentum of the 
potential flow outside, the velocity is everywhere positive and differs 
from the basic flow by a small amount in order to support the pressure 
rise. As the velocity is high, both compressibility and viscous effects 
will become important. Under this condition, the problem is considerably 

simplified by taking a uniform basic flow, namely u^^ = = 1, as 

discussed above. Accordingly, equations (20 ) reduce to 


* (1) - -P* (1) ♦ £(- (1) * i "x) 
v x (1) - -P y W * + 5 °y) 

n^ 1 ^ + p ^ = 0 

x 

< lj 




“x 


)(!) = + p^) 


(35) 


where 


H (l) = + (7 - l^u^ 

The elimination of from the first two equations gives 


( u y <1) - T x (1) ) v - E A (“y (1> - T x (1> ) 


(36) 
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After differentiating the first and the second of equations (35 ) with 
respect to x and y, respectively, an addition yields 

Ap^ = -n^ 1 ) + -ir- An^ 1 ^ 


moreover, hy definition 


H^ 1 ^ = 7M 2 P X ^ + (7 - IJmPu^ 1 ) + n^ 1 ^ 

A substitution of and AH^"^ in the energy equation leads to 


“x (l) - m 2 %x ( 1) + 2 ^( Au « (1) + 5 %x (1) ) - | An (1) - — f Z m 2 %x (1) + 
•^T + (r - D^(-V (l) - (37) 

Equations ( 36 ) and (37) form a system of linear partial differential 

equations for and v^^ with constant coefficients. The solu- 

tions, if carried out, are expressible in terms of trigonometric and 
exponential functions. According to the arguments of these functions, 
they fall into two groups: One varies slowly and is "potential-like "j 

while the other varies rapidly and, therefore, is "boundary-layer-like." 
The latter group consists of two exponentials whose arguments differ by 

a term of Hence, for large Reynolds number, the latter two expo- 

nentials will degenerate into one. It was found that this form of solu- 
tion can be obtained by solving a much simpler problem, namely, by 

assuming = 0. This assumption appears to be nothing but a method 

of approximating the solution of equations ( 36 ) and ( 37 ) in the case of 
large Reynolds number. 

If g(l) taken to be zero, then instead of equation ( 37 ) there 
is in its place 


■J 1 ) + - ^(au* 1 ) + i n^U) 


( 38 ) 
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"by elimination of pCl) and . To solve the system (36) and (38), 

two scalar functions F^^ and G ^ are introduced through the 
relations 


U^ = ^ + Gy^ 1 ^ 


> 


,(D = T Cl) 

T 


( 1 ) 


J x 


(39) 


The equivalent system^ then, is as follows: 


AG C 1 ) ■ i A AG^ 1 ^ 
A R 




-ApCl) + m^f Cl) _ af (l) _ y^Q (l) + AG Cl) 

xx x u xy r ^ u y 


According to the form of the potential solution, it will again he 
assumed that 


G (l) = g 1 (y)e Xx 

pC 1 ) = f 1 (y)e Xx - 

By the condition that G^ is finite at infinity, the solution is shown 
to he 


g-, = Co cos Xy + Cu sin Xy + — e 
^ ^ XR 


2 -cry 


( 41 ) 


where Cg, C 3 , and are constants of integration and 
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By substituting g^ from equation (4l) in the second of equations (40) 
the solution for f-j_ was found to be 


-Pi \y p-iXy 

f^_ = C 3 sin Xy - cos Xy + C^e - 1 + C^e - 1 + 


(7 - l)M^ 0^2 


R 2 + - ljxt^R 


oy 


(43) 


where 


^ - #)/(> * #) 


The velocity components for y > y,, by equation (39)? are given by 


u 


( 1 ) = 


A/C^e 


• p i x ( y - y i) ^ A x ( y - y i) 


(y - ljM 2 


E 2 + (& 


(t - 1 H R 


+ XCge 


1_ 

XR 


crCge 


-a(y-yi)l gXx 


(44) 


r(D = 


i ~h XC 5 e 


' P l X ( y - y l) + 


3-jXCge 


p i x ( y - y i) 


L 

-N 

(7 - Dm 2 , 1 

C e' 0 ( y ' yx ) > 

XR 2 + (~- - lJx^R R 

' 2 e , x 
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by redefining the constants C. It is noticed that the first group of 
exponentials varies with both x and y with a slope proportional 
to X , which, according to experimental evidence for this type of flow, 
is a small number; that is, the variation with x is relatively slow. 

On the other hand, the second exponential varies with y with deriva- 
tive proportional to g, which is of o(\/xr). For large Reynolds num- 
ber, it appears that a could be thousands of times larger than X. 
Therefore, in the case of large Reynolds number where the viscous layer 
is thin, the first group of terms changes only slightly while the second 
exponential drops to zero. Applying the boundary-layer concept, the 
potential-like terms will be taken as constant and equal to the boundary 
* value of the potential solution. Consequently, the solution can be 
written as 


( 1 ) = 

A + (k - i)oG 2 e- a ( y -M 

e Xx 


u— —— J 


> 

(1) = 

- (¥ + t)‘*" c(y " yi) 

e Xx 


(Mm) 


with 




It is seen that when cr(y - y^ becomes large, this solution joins the 
potential solution at y = 0 and the constants Cg and Cg will be 
considered as eliminated. By means of this approximate solution the 
pressure p(^) is shown to be 


P 


( 1 ) = _ 



kg 

3 


C 2 e 


-<y(y-yi) 


Xx 


(1*5) 


vhere 
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Eigenvalues X 

In the last two sections, the velocity and pressure have been cal- 
culated by entirely different methods of approximation in two different 
layers. For the incompressible layer, the nonslip conditions are sat- 
isfied on the plate, whereas in the compressible layer, the velocity 
agrees with the potential flow at infinity. The complete solution is 
then left with three undetermined constants C^, Cg, and A and with 

an arbitrary parameter X . To determine these constants, it is assumed 
that the two solutions (33) and (34) and (MO and (45) must join at the 
interface y = y^. Now, because of the simplification made in connection 

with equation (44 ) , the conditions at the interface are to be stated as 
follows for y = y 1 : 



where the subscripts c and i indicate, respectively, solutions of 
compressible and incompressible layers. The u velocity profile thus 
will have a discontinuity in slope. This could be inqproved by dropping 

the assumption = 0. However, as pressure depends only on the 

velocity and its second derivatives, an error in shear can produce only 
minor contributions and can be ignored. By substituting the solutions 
in equations (46), there results the following system of equations: 
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where U V^ 1 ^, and are defined by 


U (l) = (al) 1 / 3 


(cos TT )^ ti 1 / 2 ( cos ttj)K^ dt) + 


(sin TT ) 2 _) f 1 ^^(sin tt] )K-l dT] 
^ 0 T 


v (1) = (ax) 1 / 3 


(sin ttjJ 


(cos TTjj^ ^ J 1 Tj 1 / 2 (sin ttjJK-l dt) 


'^1 / 

t] 1 / 2 (cos TT))^ dT) 


0 


,(1) _ p o ( ° )a 


' TT1 l K 2(f \ 3/2 ) + (3^.)' 1/3 V (1) - TT li (aX)“ l/3 U (l) 


In order that linear system (47) will admit a nontrivial solution for 
Cl, c 2> and A , it is both necessary and sufficient that the determinant 



( 48 ) 
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vanish. This equation serves to determine all proper values of X. 

After , X has "been determined, the constants C-^ and C 2 corresponding 

to this X can he solved for in terms of A, namely 



where the constant A remains to he determined. 


To solve de term! nan tal equation (4-8) analytically even for this 
simplest case does not seem to he possible. The procedure from here on 
is essentially numerical# For the present purpose, the numerical solu- 
tion will he carried out, based on the approximation that tj will he 


taken mathematically as infinity and small. Accordingly, the 

integrals in and will have an infinite upper limit 

and cos tt)^ ^ 1 and sin ttj ^ % Xy^_. Furthermore, for large values 


of ru the Bessel function 


/, 3/2\ -3 A -f r, 3/2 

K llf "l S 3 1 


will he 


very small and can, therefore, he neglected. Equation (4-8) is finally 
expressed explicitly in terms of the parameters R, and a. To 

retain terms up to the order o(r"°^), 
becomes 


the determinantal equation 


Py^x 


a 1 /3ry//3 


+ X, + p 


( 0 ) 


ap = 0 


where r 



(50) 


For M = 2, R = 77k., and a = 0*332, there are two pairs of posi- 
tive and negative roots. The negative roots would make the perturbation 
infinite at negative infinity, which is contrary to assumption. Therefore, 
negative roots are not admitted. For the pair of positive roots, one 
is roughly 20 times the other which is 0 . 0 ^ 67 . For large Reynolds num- 


ber, these would correspond to X-> 


- 3/2 1/2 

»* 0 R 


and 


^2 


p 3/ V 1/4 . 
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In region 1, x is negative. The disturbance with a logarithmic rate 
of decay Xj_ will quickly disappear. The observed disturbance must 

have a decay rate equal to Xg- 

According to this solution, the dependence of the distance influ- 
enced by pressure disturbance on the Mach number and Reynolds number can 
be discussed. By the form of the first-order solution, the pressure 
decays exponentially from the point of incidence. For a given value 
of e, when -X^X = d, d being a constant, the pressure disturbance 

would have dropped to a certain fraction of its initial value. Thus, 

-X = d/Xg would serve as a measure of the distance reached by the pres- 
sure. Therefore, by varying M and R, the distance X'/L will change 

_^/k _3/k 

according to the law B ' R ' . Namely, by increasing both Mach num- 
ber and Reynolds number, the distance reached by pressure disturbances 
will decrease. This result is quantitatively different from that arrived 

at in reference 1, which is p -1 / 2 R -1 / 2 . By comparison there is an 
increase of both compressibility and viscous effects in the new result. 
This shows how gross an error can be made if the boundary-layer approxi- 
mation is applied. It is surprising also that the upstream pressure 
propagation depends only on M and R but, to the first order at least, 
is independent of the shock strength. This seems to be in agreement with 
experiment (reference 6). 

It can be shown from equations (33) and (3k) that for arbitrary 
values of A the pressure p^^Xg^jO) at the point of separation 

where Uy = 0 is proportional to p ^R - ^^ . 


FIRST-ORDER SOLUTION - DOWNSTREAM x > 0 


As in the upstream first-order solution, one must start with the 
potential solution. In order to determine the flow in region 2, the 

interaction between the incoming wave ^cp ^ ^ ^ and the shock must be 

considered. 


Potential Solution 

In equation (27 ) the incoming wave is given by 


n cp^ 1 ) = A X(x-Py) 
X 
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When this train of waves hits the shock, the shock will be slightly 
modified according to Rankine-Hugoniot conditions. The velocity poten- 
tial in region 2 with these conditions satisfied is shown to be (appen- 
dixes A and B) 



That is, to the first order, the flow in region 2 is simply a super- 
position of a step wave upon the transmitted wave. Since the conditions 
are prescribed on the shock OS, they are uniquely defined in region 2, 
terminated only by the extreme Mach line OM. 

To continue this solution into region 3; two possibilities present 
themselves. It may be either continuous or discontinuous there. If it" 
were continuous on OM, there would be a discontinuity in pressure at 
the point of incidence. In the case of inviscid flow, this, of course, 
would be admissible. But, by the condition that the viscous layer will 
join smoothly with the potential field, this would make the velocity 

u^(x,y) discontinuous at the origin x = 0. Hence, the first possi- 
bility must be discarded. If the solution is discontinuous at OM, and 

the pressure as well as the velocity u^ are continuous at the origin, 
then the discontinuity must correspond to an expansion. This, of course, 
is what has been observed. 

Assuming that the pressure returns to the value just in front of the 
shock, a simple calculation shows that the turn of the flow through a 
Prandtl-Mieyer expansion has to be 

-eP^u^Co) - 2 u(l) (°)) 

where ^u^^(O) and 2 U ^^( 0 ) are, respectively, the velocity u^) 

just before and after the shock. By means of solutions (27) and (51 ), 
the turn required is e. Wow the direction of the flow before the expan- 
sion is e (l - pA)j therefore the total inclination of the velocity 
vector at OM is 


e(2 - pA) 


A solution in region 3 subject to these conditions is found to be 
(appendix B ) : 




( 1 ) 


= i(x + py) - i(x 



k(x-py) 

e 


(52) 
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The velocity in region 3 can then he given by 


u* 1 ) = Ae x ^ x -^) 

3 v( 1 ) = 2 - pAe x ( x ~Py) 


Were the perturbation velocity to remain finite at infinity, X, would 
have to be negative in region 3; as x - Py ^ 0. 


Viscous Solutions 

As potential solution (52) in region 3 is a superposition of two 
different types of waves, namely, step waves and transmitted waves, it 
is expected that the viscous solution will also be composed of two dif- 
ferent parts, each associated with a special wave in the external field. 
Since the origin and character of the disturbances are different, they 
can best be discussed separately as the following: 

In view of the fact that the step waves are opposite in sign, they 
do not change the pressure but give rise to a uniform vertical velocity 
in the potential field. The flow in the viscous layer due to this uni- 
form deflection may vary rapidly in the y-direction but certainly not 
in the x-direction because of the constancy of pressure. The main 
effect of the step waves, first of all, then, is to produce a constant 
deflection along the edge of the viscous layer. If the velocities due 

to step waves are u g ^ and v g ^, the problem should be solved 
subject to the boundary conditions: 


u s ^) = v B (-*-) = 0 


when y = 0 


"s 


(-*-) = 0, Vg(^) = 2 at outer edge of viscous layer )> (53) 


u s (l) = 0 


when x = 0, y = 0 


It must be borne in mind that generally the flow in region 1 has 
been separated, with a considerable region of backflow in the neighbor- 
hood of the point of incidence; the resultant flow in the incompressible 

layer must be small. If Ug^^ and v g ^^ are considered as additional 
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perturbations due to the step waves, the problem can then be simplified 
by neglecting the inertia forces in the neighborhood of the plate- The 

perturbation stream function then satisfies 

A A\|r ^ = 0 (54) 

s 

in the incompressible layer OS y = y^. 


As the objective here is to demonstrate approximately the effect 
of the step waves on the viscous layer, only simple solutions will be 
considered. A simple solution that satisfies both the first and third 
of conditions (53 ) is clearly 



(55) 


for y in 0 ^ y = y^, where D is a constant. The velocities are 
thereby 


u. 


( 1 ) _ 


( 1 ) _ 



(56) 


Accordingly, the pressure, or the pressure gradients, in the incompress- 
ible layer due to the step waves is of the order R -1 and for the pres- 
ent approximation will be neglected. 

In the compressible layer where the viscous effect is not as impor- 
tant as the compressibility effect, a flow that does not associate with 
a pressure rise can either be a uniform field or the one with a vertical 
velocity varying linearly as y. According to equation (37 )> the fol- 
lowing, with the second of equations (53) satisfied, is a solution for 
y ^ y S y,) y-p being defined as the outer edge of the compressible 

layer : 
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(57) 


The fact that equation (57 ) is a solution shows the importance of both 
the compressibility and the viscosity in this layer, which allow the 

independent variation of and T^. The joining condition at 

the interface gives the constant D a value 

Thus, it is seen that the deflection of the flow due to step-waves, 
though it contributes no pressure, induces a forward velocity in the 
incompressible layer. Since it increases linearly with x, the back- 
flow should be expected to be reduced in the downstream direction. Of 
course, the possibility exists that an additional pressure might also 
enter if the inertia forces, though small, were not neglected. Never- 
theless, it can be safely stated that the pressure is always of the 
secondary importance in this case and hence the above conclusions will 
remain valid. 

On the other hand, the perturbation velocities u^"^ and v (D 

due to the transmitted wave, according to equation (52), are subject to 
boundary conditions 





To solve this problem, the viscous layer is again subdivided into incom- 
pressible and compressible layers. Since the boundary conditions and 

the differential equations are the same, and v^ ^ ^ will have 

the same forms as those given in equations ( 33 ) and (MO, with eigen- 
value X. satisfying equation (48). Now if the perturbation is required 
to vanish at positive infinity, X should be negative. But, were X 
negative, solutions (33) and (44) would be highly oscillatory, as the 
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arguments of both Bessel and the exponential functions involve a fac- 
tor \/XR. That is, the solution of the incompressible layer would 

• 3/2\ 


involve J- 


r i (|Jo 0 (0 W 


) 


and that of the compressible layer would 


have cos \/xR y and sin VxR y. For large Reynolds number, these types 
of solutions for a steady laminar flow must be rejected as impossible. 

By the same reason, complex roots are excluded. The only alternative is 
to accept the positive X. 


If X is positive, it follows that the velocity u^. (1) 

and pres- 
sure p.^^ will continue to increase until the process of linearization 

breaks down. This would imply that in the potential field the train of 
compression waves, immediately following the expansion, will grow expo- 
nentially with x. The flow direction as well as the curvature of the 
streamlines will also increase sharply. For such a flow, it is well- 
known that the Mach waves will converge and form an envelope. There- 
fore, in region 3 a shock eventually must be developed. On the other 

hand, in the viscous field, as u_^"^ + u for any constant y will 

increase with x, backflow, though slightly reduced by the step waves, 
will become stronger in the downstream direction. Now, it has been well- 
established that the laminar velocity profile with a point of inflection 
is highly unstable. As the pressure continues to grow at a Reynolds num- 
ber generally above the critical value, transition must occur after a 
critical pressure gradient is reached. The flow from there on cannot 
be theoretically studied without considering unsteady flow. 


It is therefore concluded that in the case of an incident compres- 
sion shock, laminar flow is not possible for the whole viscous layer 
and transition always occurs. This appears to be in complete agreement 
with the present available experimental observations. 


Character of Flow After Transition 

Although the flow from a certain point on is unknown, the inter- 
esting fact is, however, that the flow up to the point of transition is 
very insensitive to what happens beyond the point of transition. In 
the previous sections, the problem has been reduced to depend only on 
one constant A on which the quantitative behavior, but not the char- 
acter of the flow, depends. A quite similar conclusion has also been 
reached in reference 1 and seems to be an experimental fact (refer- 
ence 7)* In one of his experiments, Liepmann introduced an expansion 
wave immediately after the incident shock, and the observed upstream 
flow field was practically unchanged. Therefore, in order to account 
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for the observed pressure distribution, say it is sufficient just to 
take approximately the effect of the flow downstream of the point of 
transition. 

Assume that the second shock starts at the edge of the viscous 
layer and, by linear theory, has the same strength 6. Inasmuch as the 

direction of the flow on OM is constant and equal to e (2 - PA), if 

the flow after the shock TR (fig. 2) is parallel to the plate, the 
direction in front of the shock at large distance must be e. This con- 
dition gives A = p t Moreover, as, the velocity vector is turning away 
from the plate, the pressure continues to rise and will be stopped only 

by transition. In that event, the sudden increase of shear of the vis- 

cous layer will thicken the viscous layer. This thickening will make 
the flow expand and consequently a drop of pressure will ensue. After 
this , the streamlines will gradually level off and the reattachment of 
the partially turbulent layer, if not yet accomplished prior to the 
transition point, will be certain to follow. If the point of transition 
Is chosen to coincide with the point where v = 0 along the boundary 
streamline, the location is then determined by Xx-^ = log e 2. 

For values of x greater than x-^, the flow in the immediate 

neighborhood of transition would have a greater influence than that far 
downstream where the flow is more uniform. Since the flow in the tran- 
sition region is of a boundary-layer type, namely, the backflow ceases 
to be a factor, it can again be approximately represented by the inte- 
grated effects, such as the momentum and pressure. If the flows before 
and after the point of transition have the same pressure and total 
momentum, the dynamical equilibrium can then be maintained. According 
to the solution given in reference 1, the pressure distribution in the 
transition region is approximately 


!♦ 


K-> e 


XtX 


X,oX 

+ K 2 e 


where and 2 are negative constants, being functions of M 

and R and and K 2 are integration constants- By the conditions 

that at x = x-|_ the pressure and its derivative are continuous, 
and K 2 can then be determined. 

It must be emphasized again that the conditions stated in this sec- 
tion are tentative. No obvious reasons beyond the ones outlined at the 
beginning of the section can be given at this moment. 
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NUMERICAL EXAMPLES 


As a numerical example, the case of M = 2 and R = 774 is pre- 
sented. First, the determinantal equation (50) is solved numerically 
for only one Mach number 2 but at different Reynolds numbers. The 
results are shown in figure 3 • At this Reynolds number, X is 0.0467* 
The velocity u(x,y) and the surface pressure p(x,0) for this X are 
calculated. It is seen that for the shock strength of e = -1°, the 
boundary layer first separates at x se p * 1.106 (fig. 4(a)) and, subse- 
quently, backflow sets in. When the shock strength is increased to 
6 = -3°, say, the separation occurs at a much earlier station, namely 
at x ge p = -22.4. The region of backflow is proportionally wider 

(fig. 4(b)). However, if the step waves are taken into account, then 
in the case of e = -1° there would be no separation (dashed curve in 
fig. 4(a) for y^ = 3 and y£ = 9); whereas in the case of e = -3° 

the flow at the same location remains separated but the backflow has 
been very much reduced (dashed curve in fig. 4(b)). These results show 
quite the same characteristics as the experimental measurements by 
Ackeret, Feldman, and Rott (reference 8) of the velocity distribution 
over a curved plate. 

The pressure distributions over the plate for these cases are shown 
in figure 5. For the weaker shock s = -1°, surprisingly good agreement 
with experiment (reference 6) is obtained. This very fact seems to 
justify the present assumptions regarding the structure of the viscous 
layer. In the case of stronger shock, for example, e = -3°, there is, 
however, a distinct difference between theory and experiment. Theoreti- 
cally, the pressure would still decay exponentially upstream but would 
begin with a larger amplitude. Experimentally, it was found, strangely 
enough, that, over a considerable range of the upstream disturbed region, 
the pressure first decreases very slowly and then decays more or less 
exponentially. This "pressure bump" seems to be characteristic of the 
pressure distribution in the case of interaction between stronger shocks 
and the laminar boundary layer. This bumpy character in the pressure 
distribution, by all evidences, must be attributed to the nonlinear 
effect of the flow. This will be exhibited in the following section. 


APPRAISAL OF HIGHER-ORDER EFFECTS 


It has been shown that, if there is no backflow in region 1, the 
pressure disturbances will decay exponentially. When the shock strength 
increases, however, the pressure in the disturbed region becomes much 
higher and drops much more slowly than that predicted by the theory. 

This appears to be due to the fact that when backflow develops, there 
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will be an underestimate of the perturbed velocity and, consequently, 
a much lower pressure. 

To estimate the effect due to higher-order terms, a second-order 
solution is found and is given in appendix C. According to the second- 
order solution, the correction terms for wall pressure would be 

lx , 2Xx 
P-^e + Pge 

In view of the fact that there is an over estimation of the first-order 

pressure by the condition A = 3 t the correction term is expected to 
be negative for small values of x and positive far upstream. Then 

P^ > 0 and P 2 < 0. Since e^^ increases faster than x 2 and e^ x , 

when x is positive, the correction would be negatively large far down- 
stream. It has been pointed out previously that the iteration process 
fails for positive values of xj this solution shows further that it 
diverges oscillatorily. In order that the solution will behave properly, 
it will have to be carried to an odd order, such as one or three. To 
third order, the correction for wall pressure would be 


P^^ + P 2 e 2)oc + P 3 e 3Xx 

If P]_ = 0*75, Pg = -I.83, and P^ = 1, the wall-pressure distribution 

for e - -3° would resemble the curve as shown in figure 6, which does 
exhibit the same character as measured by experiments* This, of course, 
cannot be considered as conclusive evidence, but at least it shows that 
higher-order terms such as given above do have the possibility of 
accounting for the observed behavior of the wall-pressure distribution* 
It is very important that this step be carried out* 


DISCUSSION 


This study being intended as an exploratory study, the numerical 
results obtained are not expected to be exact, but only accurate enough 
to insure the conclusions* As far as the first-order solution is con- 
cerned, the entire problem depends on the determination of one param- 
eter X, which was calculated on the basis of two main approximations: 



1/2- 


¥ 



dTj 


is negligible 


(a) the contribution of the integral 
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and (t>) A,y^ is small such that terms of order and higher can 

he neglected. The subsequent solution for X satisfies this criterion. 
The question remains, however, whether Xy^ = n(l) is also a solution 

or not. If Xy ± - n(l), the integrals J ^ T ) 1 / 2 t tjj- t)3/2^qLt) 

would require a more elaborate evaluation. The possibility exists and 
is worth testing. 


As from the outset the boundary -layer approximation is rejected, 
it is of interest to examine the character of the perturbed flow in the 
light of boundary -layer theory. According to the boundary-layer approxi- 
mation, if the velocity u' = 0(1^) and 8/ck' = o(l 1 ), then the 

velocity v' = 0 fu^R and 8/8y ' = o(l~^R ). Now, from the character 


of the solution for perturbed flow in the present theory, d/dx' = 

and 8/Sy 1 = o(l~^R^^). By comparison, the derivatives for the perturbed 


flow are much larger than the corresponding derivatives of boundary-layer 
theory. Moreover, from the equation of continuity the ratio v^Vu^^ 
of the perturbed velocities is As in the vicinity of the 

wall u^ u^, v^ is much larger than the vertical component in 


the boundary layer. Since in this case 8/dx' < 8/Sy' and v^^ < u^^, 
a set of equations analogous to those for the boundary -layer flow can be 
deduced as long as backflow does not take place. Therefore, for expan- 
sion and even weak compression incident waves, a much simpler problem 
would be feasible. 


Finally, it might also be noted that, owing to the assumption that 
the flow far away from the plate is inviscid and irrotational, there is 
introduced a sharp discontinuity in higher derivatives at the demarcation 
line between the two flow fields. Because of the presence of shock, how- 
ever, a first-order discontinuity is also expected, because, by the 
assumption of frictionless flow at large distance, a smooth transition 
from small shock thickness at large distance to a larger one at the 
vicinity of the wall is precluded. The viscous -layer concept then * 
idealizes the situation by taking the shock as a discontinuity in the 
outer field but continuous in the viscous flow. The effect of this is 
exhibited in a discontinuity in slope of the streamlines. This picture 
is entirely in agreement with the observed flow patterns. 

According to the numerical example, the step waves tend to weaken 
the backflow downstream of the point of incidence but are unable to make 



ko 


) 
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the flow reattach if the separation has occurred upstream of the point 
of incidence. This may very well he due to the poor approximation 
inherent in the Stokes flow which may underestimate the rate of change 
in the y-direction. If a more exact solution is given, reattachment of 
the separated flow might he accomplished in the laminar regime. This 
point should he considered as open. 


CONCLUSIONS 


The following conclusions are drawn from an analysis hy the 
differential -equation method of the interaction of an oblique shock 
wave with a laminar boundary layer along a flat plate: 


1. The pressure perturbation decays exponentially forward of the 
point of incidence and the distance of pressu re propagatio n varies with 

different Mach and Reynolds numbers as (3 ^^R ^ where p is \Jvfi - 1, 

M is Mach number, and R is Reynolds number J and very slowly with the 
shock strength. 


2. If the shock strength is strong enough, separation of the flow 
always occurs. For given Mach and Reynolds numbers, the separation point 
depends strongly on the shock strength. 


3- The pressure at the point of separation varies with Mach and 


Reynolds numbers as 




1. In the viscous layer, la mi nar flow is not possible everywhere, 
no matter whether the incident shock is strong or weak. In the distance 
of about two or three boundary-layer thicknesses, transition would occur. 


5. The curvature of the streamlines after the shock is positive and 
the Mach waves in the potential field must coalesce to form a shock which 
approaches asymptotically the regularly reflected shock in the inviscid 
fluid. As its position depends on the point of transition, the exact 
location cannot be predicted by the present theory. 

6. From the calculated pressure distribution over the wall, it is 
definitely proved that the observed overcompression of the wall pres- 
sure is a consequence of the positive curvature of the streamlines and 
the expansion is associated with transition. 

7 . The observed "bump" in the pressure distribution in region 1 for 
strong shocks is definitely a nonlinear effect and is an indirect conse- 
quence of separation. 


Cornell University 

Ithaca, N. Y. , January 11, 1952 
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APPENDIX A 


PERTURBATION OF AN OBLIQUE SHOCK 


Let 3_ u |p p^, and P-^ denote , respectively, the parallel 

and normal velocities with respect to an oblique shock, the pressure 
and density in front of the shock, and, by replacing 1 by 2, the corre- 
sponding quantities behind the shock. These variables are related to 
each other by the Rankine-Hugoniot conditions: 


P 1 1 U J_ P 2 2 U ± 




1 U II " 2 U I 


P 1 + P 1 1 U 1 = P 2 + P 2 2 U 1 




y ii . 1 u 2 = 7 !£ 1 u 2 

r-ip 1 2 1 -L 7 - 1 p 2 221 


(Al) 


If the flow In front of shock is given, equations (Al) yield the fol- 
lowing solutions: 


Po = Pt + 


TT p i(i u l 2 ' a i 2 ) 


P 2 P 1 


7 + 1 


7-1 + 2 


> 


,i u iy 


2 U 1 l U J-^7 + 1^7 + 1 2_ u i^ 


7-1 2 a l 

1 + 


(A2) 



k2 
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where a^ is the speed of sound in front of the shock. If a ra and 
are the speed of sound and velocity at infinity, then a x can he 
expressed in terms of velocities by the relation: 





Furthermore, if oo is the angle between the velocity ( U ]/ T p) an( i the 

shock, it is easy to show that the Cartesian components before and after 
the shock are related by 



with q 2 = u x 2 + V x 2 . 

In the case of weak shock, the shock angle 0 ) can be expanded in 
terms of the deflection angle e. Then cos co and sin gd are shown 
to be 


cos o> = 7 - si + — + -- - - 

Ml ). r 3 


( r . + + x ) - 1* 

kp 1 *' |_ 8p 2 



^[i y + 1)mZ (m 2 + 1 ) - i] e 2 + 

kp 2 |_ 8p 2 
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Now since the velocity components and v^ can be expanded; 

u ± = 1 + 6u 1 ( 1 ) + . . . 

V 1 = 6V i^ + *\ (Z) + . . . 

If 03 has the expansion: 

a>-J 0 > + ,J 1 > + « 2 o.( 2 ) + . . . 

and a><°> is defined by equations (A^)^ then a straightforward reduction 
gives 

Ug — 1 4 - €U£> ^ ^ 4 - €^Ug ^ ^ 4 * « » • 

' * y 2 > evgf 1 ) + S 2 v 2 (2) + . . . (AS) 

where 
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APPENDIX B 


POTENTIAL FLOW IN REGION 2 


If the flow in region 2 is irrotational, there exists a poten- 
tial cp(x,y) and, if the shock conditions are satisfied, it must 

u 

possess the expansion: 


2 cp = x + 6 2 <ph) + e 2 + 


For the first-order solution of the shock conditions , namely , on a line 
x + py = 0, 


A (1) , [i _ _!>_/Lz A ♦ l\U> + ifi + 

2X r + 1 \ 2 m 2 / »\ (r-M)M 2 


(1) _ „ (1) 40 h - 1 i\ , (1) 


7 + i\ 2 m 2 , 


+ — VL, + 1 + 


Lp^CD 


2jD 


(7 + ljM 2 


Now on the shock 


Ul (D » Ae 2 ^ 


v 1 < 1 ) ■ -0Ae 2XjI 


If the perturbed shock angle is 


j!) = se L(x-Py) 


the solution g cp(l) must be 


cp^ 3 -) = i(x + Py) + te X ( X “ Py ) 
C P 
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b5 


Substituting cjd^ and in equation (B2), the two equations 

yield uniquely 

(B4) 


The solution in region 2 is thus a superposition of the transmitted 

waves e ^( x_ Py) a uniform step wave x + Py introduced by external 

agency. 

Since this solution will terminate on the Mach line OM, the solu- 
tion in region 3 will be found by the condition that on x - Py = 0, 

= -jU^ 1 ) (B5) 

By this condition, the Prandtl-Meyer flow requires a deflection 

8d = -ep ^^(O) - gU^^O) (B6) 

Solution (A3) gives 8d = e. Consequently, the initial direction of the 
flow in region 3; to the first order, is 

2 v (l) (0) + S3 = (2 - pA)e (B7) 
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APPENDIX C 

SECOND-ORDER SOLUTION 
Potential Field 


If the first-order potential is of the form; 


m C 1 ) = A e X(x-py) 

X 

it can he shown that the second-order equation ( 19 ) admits a solution 
which reduces , at the edge of viscous layer y = 0 , to 


u ( 2 ) _ Be ^ + Ce 2Lx 


v (2) - -SBe ’ 01 - P 


C + 


(7 + 1)mV 2 


4p‘ 


> 


2 Lx 




(Cl) 


•where B and C are arbitrary constants and the term - 


(7 + 1)mV 2 

4p 


,2 Lx 


is derived from the particular integral y e 2X.(x-|3y ). rphp 

4p 

2 

reason that hoth e and e are included in the complementary part 
of the solution is that the solution in region 2 , according to equa- 
tion (a 6 ), contains hoth solutions. 


Incompressible layer 

Assuming that,, to the second order, the density in the inner layer 
remains constant, the system of equations (15) simplifies to 


( 2 ) ~ , ( 2 ) 

A A\]r - ay At x 


= Rlt r At 


d) A * (D * (1).* (D) 


- t 


At, 


(C2 ) 
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where the second-order stream function 



is defined as 


u( 2 ) = * ( 2 ) 

y 

v (2) = -Mr (2) 
x 


By substituting + into the right-hand side of equation (C2), the 
particular integral can he foundj this is 


If 


( 2 ) _ 


sin tt) P 1 ! 

?( 2 a \) 2 / 3 ^ 0 


(cos T T| )Z ( T] ) dt) - 


cos t rj 


T ( 2 a \) 2 / 3 



(sin TTpZ(q) dif 


e 2Xx 


(C3) 


where the function Z(rf ) stands for 


m = f if3/ 2 ) 


r j 1 / 2 !^^ T)3/ 2 ^A(ti) dr) 
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with 


A(„) - i(Zffl,) 4/3 C] 


Po ( °4 


t.3/2) 


(cos Tq) r T ) 1 / 2 (cos TT1 )K 1 drj + 

Jo 3 


(sin tt]) [ T]^/ 2 (sin ttj)K^ dr] 

Jo 3 


— nK 2 (| q3/ 2 ^ (sin ttj ) r ^(cos tt))^ dT) - 
T / J 0 o 


(COS TT) ) T)^/ 2 (sin ttj)K^ dq 

J 0 3 


(C5) 


and the variables: 


f - (2SX) l/3 ^y . |£) 


(C6) 


= 2 A, /(2otX) 1 / 3 


(C7) 


The general solution, consequently, can be written as 


^ 2) = V (2) +*_ (2) 




(C8) 


f 2 ) 

where the complementary integral ^ , by the first-order solution, 

will be given by 




(C9) 
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Since the form of i|f^^(Tj) is known, and, furthermore, (0) : 

the solution i|r(2) is determined with two arbitrary constants and C^. 

With \|r(^) determined, the second -order pressure p^^(x,0) on the 
plate can be written, similarly, as 


)(^)(x,0) = P^e^* + Pge^^- + Pp e<: 


(CIO) 


vhere 


_ t^) 1 / 3 3 2/3 r n i/2 K/2„ 3 A 


P 3r 8 


)Sf/» ■> 


and P-^ and P 2 take the same form as the symbol P for the first-order 
pressure p (-*-)- 


Compres s ible Layer 


In the compressible layer, p ^ ^ = u^^ - - 1. Again, if 


r(2). 


0, equations (15) become 


( 2 >-vJ 2 >'l -Ia LS s) -vJ 2 


r A V“y 


x^^ac;^ + v^'af^ 1 ') 
^DAGyd) - uC-WyW 


mV (2) - n< 2) - ^(au ( 2) + i % ( 2 >) = -(1 + ^Irdluyd) 


(€ 11 ) 


[ 2P 2 - (l + ^)]. (l) Ty (l) - 2\(l + ’Ll! M 2 )^ 

XV?(r - ljCi 1 )) - j(''y^ 1 ^) - X v ^ 1 ^ v yy^ 1 ^ 
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Similarly, "by introducing F^ and G^ 2 ^ through the relations 


u (2) = F x ^ 2 ^ + G y ( 2 ) 

v (2) = F (2) _ q^(2 ) 


there result the following equations: 


AG. 


^ i A AG^ 2 ^ = A,(u^AG^ + v^AF^) + 


(Da r. (D _ u (l)ap.(l) 


'AGy 


M 2 ( Fxx ^ 2) + G^^) - AF^ - 2?f. a(f x (2) + Gy^ 2 ^ + i A AF^ 2 ^ 


> (C12) 


(1), (D 


U' 'V. 


-(l + 7 M 2 )vD) Uy ( 1 ) + Jj>p2 _ (1 + rM 2) 

2X^1 + — - " - 1 ^u( 1 ^ + \(y - 1)M 2 vD) 2 _ 1 V y(D 2 - i v^ 1 Vyy( 1 ^ 


(2) (2) 

The solutions F v ' and G v ' can he written in an analogous 


manner. That is: 


p(2 ) _ F (2 ) + ji (2 ) 

C P 


0 < 2) = Gc ( 2 > + G < 2 > 


(013) 


(2) (2) 

Here Fq and Gq are known from the first-order solution and 

( 2 ) ( 2 ) 

Fp v and Gp are given as follows: 


F p (2) = ( K 1 + K 2 e_0y 

G p ^ = K^e _0y+2Xx 


+ K 3 e 


-2oy\ 2Xx 


(C1*0 
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where the constants K are defined hy 

1 + (p 2 /^ 2 ) + (7 - D* 2 ^ 2 + (7 - I)** + 2 + 2P X 4 a2 + 

1 8x(p x 2 - p 2 2 ) + 8 J**. 

3 

x - (p 2 /p x 2 ) [2P 2 - (7 - i^Pj 2 + 2 + (7 - i)t^ 

8xp ? 2 . Jyift 

2 u jT 


1 - (P/Px) 


J l * 2 + (xp x - a) 2 ] 

’ % ' 47 |$)L r( “ 1 '" ) L f 55 — ~J X 

(^ti 2 - - k " 1 ♦ p 1 ( t(k + ir 1 )x- 2 j 

8x 2 (i - 2XR -1 ) + (XP X - of + 2(1 - ItXR" 1 ) - (XR)* 1 (XP 1 - of 


Pj(i + 7&) - ^ jit + 2(7 - 1)*^ + Pj 2 ^ + (7 - 1)»^ - P 2 ]} • 


(fff- - R -1 j ^Xp-^Pj 2 + ^ + 2J? + 2^ + a jl + 71^ + Pj 2 [(7 - 1)mP - P 2 ]^ 


CgA 


(xp 1 + of - kx\ z 


1 + (P/Pl) 1 

WxPx + 

[-- + (XR)” 1 ] 

L r 1 

[U 2 + (XP 1 - a) 2 ]] 

J-(>V • A 

k ^-I-¥ ( k+r '^ 

4 * 1 4 ) 

l 

8 x 2 (l - + 2 (xp 1 + <r) 2 (l - 

f) - (xRrl ( xp x + a ) U 


P x (l + 7 ^) - ^5 + (7 - 1)*^ - P^] + k + 2 (7 - + 

^ - R- 1 ^-Xp 1 (p i 2 + o 2 + 2M 2 + 2) + ||l + 7^ + P x 2 ^[(r - DM 2 - 



X 


rH|W 
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and the constants K and (3^ stand for 


R + 7 - lW\ 




8 7^X 

3 R 


(C15) 


The second-order solution now involves six arbitrary constants and,* by 
the conditions at the interface,* all but one can be determined. 
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(a) € = -1°J M = 2; R = 77^. 


Figure 4.- Velocity distribution. 
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(b) e = -3°; M = 2j R = 774 . 
Figure 4.- Concluded. 
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